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Abstract 


Dust is one of the natural hazards that occur mainly in arid and semi-arid areas of the world such 
as Iran. In recent years, dust storms have been one of the most important air pollution crises in 
Hamadan province. Identifying the dust and sand sources is the first step in controling this 
phenomenon. Thus, this study tried to identify dust sources of Hamadan province using Aerosol 
Optical Depth (AOD) from the MODIS satellite images and numerical modeling to predict 
airflow from 2008 to the end of 2018. In the next step, using HYSPLIT model, the path of 
particles’ moves and their source were determined at three levels of 500, 1000 and 1500 meters 
for 48 hours before the storm for external sources and 18 hours for internal sources. The 
meteorological data of the model was obtained from a 0.5-degree GDAS. Finally, its movement 
from source to Hamadan urban area was monitored using GIS and Spatial Analysis Tools. The 
results showed that the origin of the dust phenomena in this city are three sources including 
internal parts the province (Famenin, Malayer and Razan county), outside the province (such as 
Kurdestan, Markazi, and Khuzestan province) and out of Iran (Western parts of Iraq, Syria, Jordan 
and northern Saudi Arabia). So, it is necessary to monitor the spatial patterns of dust sources in 
the western regions of Iran for different years. The combination of remote sensing data, Hysplit 
model, and field based data could well show dust transport corridors. 
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مقاله پژوهشی 


مسیریابی حرکت گردوغبار بر اساس غلظت ذرات معلق و مدل هواشناسی در استان همدان 
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حکیده 


گردوغبار یکی از مخاطرات طبیعی است که عمدتاً در مناطق خشک و نیمه‌حشک جهان ازجمله ایران رخ 
می‌دهد. در سال‌های Gol‏ طوفان‌های گردوغبار به یکی از مهم‌ترین چالش‌های آلودگی هوای استان 
همدان تبدیل شده است و شناسایی کانون‌های تولیدکننده این پدیده اولین گام در مدیریت و JES‏ آن به 
شمار می‌رود. ؛ ازاین‌رو در این مقاله به شناسایی کانون‌های تولید گردوغبار در استان همدان با استفاده از 
شاخحص عمق Sal‏ اخذشده از تصاویر سنجنده مادیس و همچنین» مدل‌سازی Gade‏ پیش‌بینی جریان 
هوا در بازه زمانی سال ۲۰۰۸ تا پایان سال ۲۰۱۸ پرداخته شده است. مسیر باد در £A‏ ساعت قبل از طوفان 
برای Lin‏ خارجی و ۱۲ ساعت برای Lite‏ داخلی تعیین شد. داده‌های هواشناسی مدل نیز از سامانه 
یکپارچه‌سازی داده‌های جهانی دریافت شد. درنهایت با استفاده از سیستم اطلاعات مکانی و ابزارهای 
تحلیل مکانی» نمايش حرکت گردوغبارها از Lax‏ شکل‌گیری تا fel‏ محدوده استان همدان ردیابی شد. 
نتایج نشان داد که Lex‏ پدیده‌های گردوغبار ورودی به این استان دارای سه منشا داخل استان 
(شهرستان‌های فامنین ملایر و رزن» خارج از استان (استان‌های کردستان. خوزستان و مرکزی) و خارج از 
کشور (غرب GLE‏ سوریه. اردن و شمال عربستان) است که بیانگر لزوم مطالعه پایش مکانی 
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گردوغبارهای این منطقه است. درنهایت ترکیب روش‌های سنجش زدوری و مدل های‌اسپلیت به همراه 
اطلاعات زمینی گردوغبار می‌تواند دالان‌های تحت تأثیر گردوغبار را به‌حوبی نشان Jas‏ 


کلیدواژه‌ها: گر NET le ss‏ عمق اپتیکی آئروسل» همدال. 


۱- مقدمه 

تحقیقات نشان داده است که قدمت پدیده گردوغبار به بیش از ۷۰ میلیون سال می‌رسد OLE)‏ و «OUS‏ 
۶ سازمان جهانی هواشناسی نیز پدیده گردوغبار را کاهش سطح دید به زیر ۱۰۰۰ متر در اثر غبار برخاسته از 
بادهای شدید تعریف می‌کند (تقوی" و همکاران, ۲۰۱۷) که دارای تأثیرات منفی بسیار بزرگ اجتماعی-اقتصادی» 
بهداشتی و محیطزیستی است (ژانگ " و هکاران, ۲۰۱۵)؛ به‌طوری‌که گردوغبار منجر به ایجاد اختلال در توازن 
تشعشعات. دمای هوا؛ مشخصات اپتیکی» طول عمر ابرها (ونگ و Lus‏ ۲۰۰۵ فرآیند بارش ژئوشیمی خاک 
نواحی تحت SL‏ (مندز* و همکاران, ۲۰۰۷)» چرخه‌های بیوشیمیایی محیط‌های آبی» فرآیندهای هیدرولوژیکی و 
ویژگی‌های زیبایی‌شناختی می‌شوند. ذرات گردوغبار علاوه بر تأثیرگذاری بر شرایط جوی و تولیدات کشاورزی» بر 
سامت انسان و اکوسیستم‌ها نیز اثرگذار هستند (هانبرگ و نیکول'» ۲۰۱۲) به‌طوری‌که در بعضی از موارد ذرات 
گردوغبار موجب انتقال باکتری‌هاء ویروس‌ها و قارج‌های بیماری‌زا نیز می‌گردد (جانسون" و همکاران ۲۰۱۱). 
همچنین با تأثیر بر محصولات کشاورزی» خاک علف‌های هرز, OUI‏ و بیماری‌ها به‌طورکلی کش‌اورزی و امنت 
غذایی را تحت تأثیر قرار داده است که خود منجر به مهاجرت به شهرها و ایجاد مشکلات بعدی شده است (ممتاز و 
همکاران. ۲۰۱۹). به‌طورکلی» اکوسیستم‌های خشکی و سلامت جوامع تحت تأثیر گردوغبار قرار دارند (فاربردر و 
همکاران, ۲۰۱۷؛ حیدریان" و همکاران» ۲۰۱۸). پدیده گردوغبار بر اساس شدت و اثرگذاری بر سطح دید افقی به 
سه دستۀ طوفان گردوغبار " (بادهای شدید و کاهش دید افقی از صفر به ۱۰۲ متر)» گردوغبار وزشی "" (بادهای 
ملایم‌تر و کاهش دید افقی تا "۱۰-*۱۰) و گردوغبار معلق "' (کمترین شدت و کاهش دید افقی به کمتر از ۱۰۴ متر) 
تقسیم‌بندی می‌شود (کنه " و همکاران (ASA‏ با وجود اینکه گردوغبار یک پدیده طبیعی و شایع در مناطق SES‏ و 
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سال یازدهم مسیریابی حرکت گردوغبار بر اساس غلظت Ww ne‏ 


نیمه‌خحشک است CSS)‏ و همکاران. ۲۰۱۵) اما در برخاستگاه نشان‌دهنده یک فاجعه طبیعی است (تقوی و OUS‏ 
۷ و ایجاد کانون‌های تولید گردوغبار نشان‌دهنده وخامت و وضیعت قهقرایی محیط‌زیست در آن محدوده است. 

در سال‌های اخیر فعالیت‌های انسانی مانند جنگ مشکلات اقتصادی, مدیریت ناصحیح منابع طبیعی» تالاب‌ها و 
دریاچه‌ها همانند تغییرات اقلیمی بر روی پوشش زمین تأثیرات گسترده‌ای داشسته و باعث تشدید خشکسالی‌ها و 
کاهش ذخیره آب و رطوبت در برخی از کشورها بخصوص ايران و همسایگان غربی آن مانند سوریه و عراق شده 
است (نوری" و همکاران ۲۰۱۹). ایران علاوه بر اینکه در نواحی خشک و نیمه‌عشک جهان قرار گرفته که خود 
موجب ایجاد مناطق بیابانی شده است» در مجاورت آن نیز منابع اصلی ریزگردها مانند پیابان‌های عراق» عربستان و 
اردن قرار دارند (باغبانان " و همکاران. ۲۰۲۰). در یک دهه گذشته گردوغبار به یک پدیده غالب و شايع تبدیل شده 
(باغبانان و همکاران ۲۰۲۰) و استان همدان یز در ارتباط با این پدیده وضعیت بحرانی داشته است. استان همدان در 
سال‌های اخیر کاهش در میانگین سالانه رطوبت. میانگین عمق برف و میانگین بارش سالیانه و افزایش در میانگین 
یر مالا و تاد ووزهای همراه با le‏ را ته کرده ات QUA UA a’ ead‏ بااین وجرد همه انارو 
نشانه‌های این پدیده پیش‌رونده» پژوهش‌های کافی در این زمینه بخصوص در مقیاس محلی انجام نشده است و این 
تغییرات اقلیمیسمحبط‌زیستی آثار ویرانگر خود را روی مردم این نواحی گذاشته است؛ ازاین‌رو شناسایی بر خاستگاه 
گردوغبار برای مقابله با این پدیده حائز اهمیت است و منجر به اتخاذ رویکردهای مدیریتی می‌شود (اسماعیلی "و 
همکاران (Qe‏ 

مطالعات مشابهی نیز در خصوص بررسی وضعیت و منشأً یابی طوفان‌های گردوغباری در plo‏ استان‌های ایران 
صورت گرفته است. ازجمله فاربایی و همکاران (۱۳۹۸) با استفاده از سنجه غبار آلودگی و داده‌های طیفی سنجنده 
Lal ob olds uasa‏ گردوغبار هو ترا غرټ و وب غرب آیران ded te by‏ کار ان Or)‏ ند جر 
منطقه جازموربان طوفان‌های گردوغباری را به لحاظ مکانی و زمانی پایش و بررسی کردند. دانیالی و همکاران 
(۱۳۹۷) نیز جهت تحلیل مکانی گردوغبار استان خوزستان از تصاویر مودیس» شاخص عمق اپتیکی آثروسل و 
مدل های‌اسپلیت استفاده کردند. پارمرادی و همکاران (۱۳۹۸) مسیرهای ورود طوفان‌های گردوغباری در نیمه شرقی 
ین را با استفاده از مدل های‌اسپلیت ردیابی کردند. جبالی و همکاران (۱۳۹۹) تغییرپذیری مکانی دید افقی leta‏ 
درگیر با طوفان‌های گردوغبار در استان یزد را با سنجنده مودیس و الگوریتم دیپ‌بلو ارزیابی کردند. رایگانی و 
خیراندیش OY)‏ با بهره‌گیری از سری زمانی داده‌های ماهواره‌ای به اعتبارسنجی کانون‌های گردوغبار در استان البرز 
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خحصوصیات فیزیکی- شیمیایی گردوغبار در منبع تولید وابسته به خصوصیات سطح زمین است (ونگ " و همکاران, 
۵ و این رویکرد منجر به درک مکانیسم‌های تولید و انتشار گردوغبار و انتخاب بهترین استراتژی مدیریتی منطقه 
می‌شود WS)‏ و همکاران. ۲۰۱۵). بر اساس مطالعات صورت گرفته در سال‌های اخیر» روش‌های مختلفی جهت 
pl tis‏ کانون‌های گردوغبار ازجمله تحلیل داده‌های هواشناسی (زو" همکاران, ۲۰۲۰) فن‌آوری RS-GIS‏ 
(بلورانی * و همکاران, ۲۰۲۰ اطلاعات زمین‌شناسی Paus)‏ ۲۰۱۳) و مدل‌سازی (اشرفی" و همکاران, ۲۰۱۶) استفاده 
شده است. بااین حال در پژوهش‌های گذشته این نکته وجود دارد که در بیشتر آن‌ها منبع گردوغبار در سطح منطقه‌ای 
و جهانی بررسی و شناسایی شده است» درصورتی که این رویکرد نمی‌تواند پاسخگوی نیازهای مدیریتی در سطح 
محلی باشد (حیدریان و همکاران» ۲۰۱۸). در این راستاء با توجه به نیاز مبرم بخش‌های مدیریتی به درک شرایط و 
وضعیت بحرانی محدوده موردمطالعه. برای اولین بار با یک رویکرد جامع در سناریوهای مختلف» شناسایی کانون‌ها 
و مسیرهای تولید گردغبار در داخل استان همدان موردبررسی قرار گرفت. 

سنجنده مادیس با دارا بودن ۱١‏ باند حرارتی برای پایش گردوغبار در سراسر جهان کاربردی و قابل‌استفاده است. 
هر تصویر رده آثروسل سنجنده مادیس حاوی طیفی از شاخص‌های عمق اپتیکی " در طول‌موج‌های مختلف بوده که 
عموماً توسط الگوریتم‌های دیپ بلو“ هدف‌های تیره" و یا ترکیب آن‌ها تهیه گردیده است. باندهای ۱-۷ برای مطالعه 
آثروسل‌ها تعبیه‌شده و باندهای ۰۱ ۳ و ۷ در بازیابی اطلاعات گردوغبار پرکاربردتر هستند. شاخحص عمق اپتیکی 
آئروسل نیز یک پارامتر کمی بدون بعد جهت نشان دادن تراکم و تمرکز ذرات موجود و معیاری از ممانعت ذرات 
معلق گردوغبار در برابر رسیدن نور خورشید به زمین است (رجبلو و همکاران» ۱۳۹۷). این شاخص در طول یک 
مسیر عمودی بر روی تصاویر مادیس محاسبه می‌شود و ویژگی‌های فیزیکی و شیمیایی ذرات» اندازه و روشنایی 
قابل‌اندازه گیری اسست. مقدار عمق اپتیکی می‌تواند با تراکم تعداد آثروسل‌ها و ویژگی‌های of‏ ذرات متفاوت باشد 
(ارجمند و همکاران, ۱۳۹۷). به‌طورکلی» سنجنده مادیس دارای چهار محصول برای تعیین AOD‏ است. 

پروداکت Aerosol Optical Depth (Ikm)‏ اطلاعات روزانه چگالی گردوغبار را برای نواحی آپی و Sse‏ با 
توان KS‏ یک کیلومتری برای محصول و توان SSB‏ ۱۰ کیلومتری برای الگوریتم ارائه‌شده در اختیار کاربر قرار 


می‌دهد که از سال ۲۰۰۰ برای ماهواره ترا وسال ۲۰۱۲ برای آکوا در دسترس است. 
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پروداکت Aerosol Optical Depth (3km)‏ از سال ۲۰۰۲ برای ماهواره آکوا و سال ۰ برای ماهواره ترا در 
دسترس است. توان تفکیک ۳ کیلومتری الگوریتم ارائه‌شده در این محصول امکان مطالعه گردوغبارها را در مناطق 
شهری به‌صورت بسیار مطلوب و کاربردی ارائه می‌دهد. 

پروداکت Merged DI/DB AOD‏ با ترکیب دو الگوریتم DT.‏ و DB.‏ امکان مطالعه بهتری از تغییرات AOD‏ را 
در مقیاس جهانی و در هر دو محیط Se‏ و آبی فراهم می‌آورد که از سال ۲۰۰۰ برای ماهواره ترا و آکوا در 
دسترس است. این پروداکت با توان SSH‏ دو کیلومتری تصاویر و ۱۰ کیلومتری الگوریتم برای نواحی حشکی؛ 
اقیانوسی و پوشش گیاهی مورداستفاده قرار می‌گیرند. 

پروداکت Deep Blue AOD‏ از سال ۲۰۰۰ در دسترس اما تنها برای نواحی خشکی قابل استفاده است. مقدار 
آثروسل توسط الگوریتم DB‏ برآورد شده که مختص نواحی خحشکی و بیابانی است. توان تفکیک این الگوریتم ۱۰ 
کیلومتری و توان تفکیک پروداکت‌ها دو کیلومتری است. 

از طرفی دیگرء ازجمله مهم‌ترین مدل‌های شناسایی مناطق خاستگاه و پیش‌بینی گردوغبار می‌توان به مدل‌های 
DREAM NAPPS COAMPS .CALIPSO WRF HYSPLIT RegCM 5‏ اشاره کرد. مدل پخش 
اتمسفری‌های اسپیلت! جهت شبیه‌سازی خطوط سیر گردوغبار و شناسایی برحاستگاه‌ها و کانال‌های گردوغبار به کار 
گرفته می‌شسود که از مدل‌های معروف جهت محاسبات انتقال جرم در اتمسفر است و کاربرد گسترده‌ای جهت 
مدل‌سازی‌های مربوط به نحوه پخش و پراکنش بسته‌های هواو محاسبات غلظت دارد. از ویژگی‌های این مدلء 
قابلیت استفاده از داده‌های هواشناسی جهانی. مدل‌سازی انتقال آلاینده‌ها در فواصل بسیار زیاد (بیش از ۱۰۰ کیلومتر) 


۲- مواد و روش‌ها 

۱-۲- منطقه مطالعاتی 

استان همدان در گستره‌ای به مساحت ۱۹۶۱ کیلومترمربع» منطقه‌ای کوهستانی و همراه با رودخانه‌های متعددی 
است که در غرب ایران قرار گرفته است. آثار ch JL sie‏ چندساله و تشنگی زمین در استان همدان در سال 
۸ آغاز شده و مناطق شمال غربی تقریباً جزء مناطقی هستند که با خشک‌سالی مواجه است. وضعیت AS‏ هوای 
استان همدان توسط پنج دستگاه سنجش کیفیت هوای اداره محیط‌زیست در ایستگاه‌های همدان, نهاوند. Ue‏ 
کبودرآهنگ و اسداآباد موردسنجش قرار می گیرد. اطلاعات هواشناسی از ایستگاه‌های سینویتیک همدان, ملاین 
تویسرکان» نهاوند. اسدآباد فامنین و رزن دریافت و برای تکمیل اطلاعات هواشناسی از ایستگاه‌های مجاور و خارج 


از استان نیز استفاده می‌شود. همچنین» ce‏ تحلیل پدیده گردوغبار در سطح استان» محدوده‌ای فراتر از مرزهای 


1 HYSPLIT 


Lil a> ۱۹۰‏ و مخاطرات محیطی شمارة اول 


v‏ استان همدان (محدوده 1٥0‏ کیلومتری) به‌عنوان محدوده اثرگذار در نظر گرفته شد که بخش‌هایی از استان‌های 


کردستان. کرمانشاه لرستان» مرکزی» زنجان و قزوین را شامل می‌شود (شکل Q‏ 


ایستگاه بابش آلودگی pu E‏ 
u-— n‏ 
مهد وده مطالعانی E‏ 
۶ 80 40 20 0 


بمپبهبپ] 


شکل ۱- منطقه موردمطالعه؛ موقعیت ایستگاه‌های سینوپتیک و سنجش ذرات Glas‏ در محدوده استان همدان 


TT‏ روش پژوهش 

دریافت cole bl‏ مورد نیاز از ایستگاه‌های سینویتیک و سازمان محیط‌زیست 

در ابتدا روزهای همراه با گردوغبار از ایستگاه‌های هواشناسی و سپس داده‌های غلظت ساعتی ذرات معلق 
Pm10)‏ و (Pm2.5‏ در ایستگاه‌های سنجش و پایش کیفیت هوای سازمان حفاظت محیط‌زیست برای ۱۰ سال اخیر 
(۲۰۰۸-۲۰۱۸) در محدوده استان همدان دریافت شد. موقعیت و مشخصات ایستگاه‌های سینوپتیک و سنجش ذرات 


معلق در حدول )۱( نشان داده شده است. 


حدول ۱- موقعیت ایستگاه‌های سینوپتیک و سنجش ذرات معلق در استان همدان و سال شروع ثبت اطلاعات 


P 
sl 


ایستگاه سینوپتیک-سال شروع ثبت داده- موقعیت جغرافیایی یستگاه سنجش ذرات معلق- سال شروع ثبت داده- موقعیت 


ایستگاه جغرافیایی ایستگاه 
نام ایستگاه سال طول عرض نام ایستگاه سال طول عرض 
همدان YEA Nor You‏ همدان No Y‏ ۳۹/۸ 
نهاوند AH Nt Yev‏ نهاوند LAY YAY‏ ۳/۹ 
EMAY 1۳ m AÛ N10 You m"‏ ۳/۹ 


ایستگاه سینوپتیک-سال شروع ثبت داده- موقعیت جغرافیایی ایستگاه سنجش ذرات معلق- سال شروع ثبت داده- موقعیت 
ایستگاه جغرافیایی ایستگاه 
نام ایستگاه سال طول عرض نام ایستگاه سال طول عرض 
تویسرکان LAY 1۳ Kals | ۰ NEY Yoo‏ 0/11 
اسدآباد YEN N۱۱ YA‏ اسدآباد ۳۰۵ Y£NA ÉN'NY‏ 
فامنین ۱ NAA‏ ۳0/۱۱ 
ون 1۱ \+/£4 Yoro‏ 


تعیین سناریوهای موردبررسی جهت شناسایی کانون‌ها و مدل‌سازی جهت > CS‏ گردوغبار 

سناریو اول- شامل روزهایی از ۱۰ سال اخیر که مقدار هر دو شاخص ۳۱۷۲10 و PM2.5‏ بالاتر از استاندارد است 
(استاندارد PMIO‏ برابر ۱۵۰ و استاندارد ۳۷12.5 برابر ۳۵ است) و هم‌زمان توسط اداره هواشناسی نیز به‌عنوان روز 
آلوده (قابلیت دید زیر ۰ متر) ثبت شده است. سناریو دوم- شامل روزهایی از ۰ سال اخیر است که در آن مقدار 
هر دو شاخص ۳10 و ۳۸2.5 به‌طور هم‌زمان بالاتر از استاندارد هستند. سناریو سوم- این سناریو شامل روزهایی 
از ۱۰ سال اخیر است که به‌صورت پیوسته (بیش از دو Gay‏ دارای مقادیر بالاتر از استاندارد در هر دو شاخص 
0 و Pm2.5‏ هستند. 

استفاده از تصاویر ماهواره‌ای 

در مرحله بعد اقدام به دانلود تصاویر سنجنده مادیس از ماهواره‌های ترا و UST‏ مربوط به روزهای آلوده در سه 
سناریو فوق گردید. جهت تعیین چگالی ذرات معلق و همچنین تعیین کانون‌های احتمالی گردوغبار بر روی تصاویر 
مادیس از پروداکت Aerosol Optical Depth (3km)‏ استفاده و برای هر روز دو تصویر دانلود شد. استخراج اطلاعات 
مربوط به شاخحص AOD‏ با استفاده از نرم‌افزار 5.3 ENVI‏ و افزونه MCTK‏ انجام شد. از الگوریتم Optical Depth‏ 
Land and Ocean‏ جهت بازسازی غلظت ذرات Glee‏ استفاده و خروجی به فرمت dat‏ تهیه شد. سپس چگالی 
ذرات Glee‏ در هر تصوير مشخص شد. مقدار شاخحص AOD‏ در محدوده ۰/۱-۰/۵ است؛ ازاین‌رو محدوده ۰/۲- 
۱ به‌عنوان شاخص آسمان صاف و کاملاً آبی و بدون گردوغبار و شفاف و با حداکثر قابلیت دید در نظر گرفته 
می‌شود. با افزایش شاخص از ۰/۱ به سمت ۰/۵ و ۱ میزان تراکم آثروسل نیز افزایش یافته و قابلیت دید PAT‏ و 
عمق اپتیکی افزایش می‌یابد. فراوانی عمق اپتیکی بالاتر از ۰/۵ نیز به‌عنوان معیاری از وقوع طوفان گردوغباری در JE‏ 
گرفته می‌شود CLAU)‏ ۲۰۱۹ بنابراین» آستانه ۰/۳ به بالا Sly‏ آسمانی کدر و تار با حجم متراکمی از آثروسل‌ها در 


جو و حداقل دید است؛ ازاین‌رو palis‏ بالاتر از ۰/۳ بر روی تصاویر به‌عنوان منشأ و کانون گردوغبار در نظر گرفته 


1 NASA 


MY‏ جغرافیا و مخاطرات محیطی شمارة اول 


شدند. در شکل (۲) تصویری از چگونگی بازسازی میزان غلظت ذرات با استفاده از تصاویر ماهواره‌ای مادیس بر 


روی col‏ ناسا نشان داده شده pga) Eal‏ مربوط به NY‏ ژوئن ۶ است». 


۷۹۰۶ WORLDVIEW 


OPTICAL DEPTH vg 


Aprono Optical Depth Iim (Jam, Land 
and Ocean) 
Agua / MONS 


400 


Start Comparison 


شکل ۲- بازسازی گردوغبار با استفاده از سنجنده مادیس؛ با حر کت از طیف رنگی زرد به سمت قرمز شاخص 


AOD‏ افزایش می‌یابد 


مدل‌سازی به روش های‌اسپلیت 

به‌منظور بررسی خط سیر ذرات گردوغبار در مدت‌زمان وقوع Ob yb‏ گردوغباری در ۱۰ سال اخیر نیز از مدل 
های‌اسپلیت استفاده شد. مدل های‌اسپلیت مدلی دوگانه برای شناسایی مسیرهای گردوغبار پراکندگی و ته‌نشینی آن با 
استفاده از رویکردهای PUFF‏ و ذرات و با استفاده از حداقل پارامترهای هواشناسی است. این مدل یک سیستم مبتنی 
بر وب است که توسط اداره ملی اقیانوسی و جوی آمریکا! در سال ۱۹۹۷ راه‌اندازی گردید. شیوه محاسباتی این مدل 
ترکیبی میان دیدگاه‌های اویلری " و لاگرانژی " است (شن * و همکاران ۲۰۰۹). مناسب‌ترین رویکرد برای مسیرهای 
گردوغبار رهیافت لاگرانژینی است. به‌طوری که هر ذره در زمان‌های مختلف موقعیت‌های متفاوتی خواهد داشت؛ 
بنابراین رهیافت لاگرانژینی این OLS‏ را pal‏ می‌کند که مدل‌های پایه لاگرانژی پس از انتخاب ذره در موقعیت 


مکانی و زمانی حاصء > ES‏ آن را در گام‌های زمانی دلخواه به سمت پیش‌بینی آینده روش * (پیشرو) و يا به سمت 


] NOAA 

2 Eulerian 

3 Lagrangian 
4 Shen 

5 Forward 


زمان‌های گذشته (روش پسگرد) ردیابی نمایند (اسکادرو " و همکاران» ۲۰۰7). برای ردیابی مسیر ورود گردوغبار در 
روزهای همراه با این پدیده و همچنین شناسایی کانون‌های گردوغباری. با مراجعه به تارنمای NOAA‏ داده‌های 
روزانه eg‏ شده توسط این سازمان با فرمت BIN‏ دریافت و با استفاده از مدل های‌اسپلیت مسیر حرکت ذرات و 
Lisa‏ آن‌ها در سه ارتفاع مختلف از سطح ایستگاه موردنظر (۰۱۰ ۲۰۰ و ۵۰0 متری) مورد ردیابی قرار گرفت. داده‌های 
هواشسناسی مدل نیز از سامانه یکپارچه‌سازی داده‌های جهانی " و نیم درجه به‌عنوان ورودی مدل انتخاب شدند. با 
روش پسگرد. هم‌زمان با آغاز اولین گردوغبار در منطقه موردمطالعه, مسیر باد در EA‏ ساعت قبل از طوفان برای Lise‏ 
M DE‏ قل او رفن بان طاشن کی E‏ بای نمی ETE NET‏ 
خروجی آن با نقشه‌های همدیدی و تصاویر ماهواره‌ای مقایسه شد. خروجی مدل Ua‏ شکل‌گیری و نحوه‌ی انتشار و 
CS >‏ طوفان‌ها را با دقت بسیار زیادی بازیابی می‌کند. 


۳- نتایج و بحث 


۱-۳- تعیین روزهای آلوده و تصاویر ماهواره‌ای در هر سه سناریو 
بر اساس سناریو Syl‏ در ۵ ایستگاه همدان, نهاوند. ملایر» تویسرکان و اسدآباد ۳۰۷ روز دارای مقدار بالاتر از 
استاندارد در هر دو شاخحص PMIO‏ و PM2.5‏ بودند (ایستگاه همدان ۱٩‏ روز نهاوند 17 روز ملایر ۸۷ روز و 
کبودرآهنگ ۱۳۵ روز). بر اساس گزارش اداره هواشناسی نیز تعداد رخدادهای ریزگرد گزارش‌شده با قابلیت قدرت 
دید کمتر از ۱۰۰۰ متر 1۲۰ روز شناسایی شد (شکل ۳. تعداد روزهای آلوده به ترتیب برای سال‌های ۲۰۰۸ ۲۰۰۹ 
۰ ۲۰۱۱۱ ۲۰۱۱۲ ۲۰۱۱۳ ۲۰۱۸ ۲۰۱۵ ۲۰۱۳ ۲۰۱۷ و ۲۰۱۸ به ترتیب NY ۳۳ ۸۵ TV Y «07 00 EA‏ 
۷ و VV‏ روز ثبت شده است. با مطابقت اطلاعات هواشناسی و محیط‌زیست تنها ۲۶ روز در هر دو ارگان به‌عنوان 
روز آلوده ثبت شده است. در سناریو دوم. ۳۰۷ روز گردوغباری در ۱۰ سال اخیر دارای ذرات UL Glee‏ از استاندارد 
بوده و توسط دستگاه‌های سنجش کیفیت هوای اداره محیط‌زیست ثبت شده است (ایستگاه همدان VA‏ روز نهاوند 
7 روز ملایر ۸۷ روز کبودرآهنگ ۵ روز). در سناریو سوم ۱ روز آلوده ثبت شده است که در ایستگاه همدان 
٤‏ رون در ایستگاه نهاوند £4 روز (به‌صورت ۱٩‏ دوره چندروزه) ایستگاه ملایر VÀ‏ روز (به‌صورت 1 دوره چندروزه) 
و ایستگاه کبودرآهنگ ۹٩‏ روز (به‌صورت ۱۲ دوره چندروزه) مشاهده شده است. جهت تهیه نقشه رقومی شدت 
غلظت ذرات در برخاستگاه‌های استان و محدوده اثرگذان درمجموع برای سناریو اول EA‏ تصویر برای سناریو دوم 
۰ تصوير و برای سناریو سوم ۳۲۲ تصوير از ماهواره‌های ترا و آکوا دانلود و پردازش شد. pl‏ تصاویر به دلیل 

پوشش ابر و يا عدم وجود ذرات Glee‏ بر روی تصاویر حذف شدند. 
Backward‏ 1 


2 Escudero 
3 GDAS: Global Data Administration System 


4٤‏ جغرافیا و مخاطرات محیطی شمارة اول 


۲-۳- غلظت ذرات glee‏ شناسایی برخاستگاه و مدل‌سازی > CS‏ طوفان‌های گردوغبار با تصاویر مادیس 

۱۰ بر روی تصاویر پردازش شده در بازه زمانی‎ (AOD مربوط به بازسازی غلظت ذرات معلق (شاخحص‎ als 
داد که چگالی ذرات در محدودۀ ۰/۹۸- ۰/۲ متغیر است. تراکم آثروسل زیر ۰/۳ نشان‌دهنده آسمانی صاف‎ OLE ساله‎ 
گردوغبار و تراکم بالاتر از ۰/۵ بیانگر وقوع گردوغبارء کاهش دید و افزایش عمق اپتیکی است. برآیند‎ BE و شفاف و‎ 
ذرات برای تمام‌روزهای گردوغباری در هر سه سناریو به‌عنوان مسیر حرکت طوفان‌های گردوغبار‎ OS > مسیر‎ 
تعیین شد. نتایج تحلیل تصاویر ماهواره‌ای در سناریوهای مختلف و مدل‌سازی روند حرکت ذرات گردوغبار با مدل‎ 
های‌اسپلیت نشان داد که به‌طورکلی تمامی مسیرهای گردوغبار استان همدان از سه الگوی داخل استان, حار استان‎ 
طوفان‌های‎ ES > و خارج از کشور تبعیت می‌کنند. -خطوط رنگی (آبی» قرمز و سبز) مسیرهای پیش‌بینی‌شده‎ 
گردوغبار را در تاریخ‌های مشخصی نشان می‌دهند؛ به‌طوری‌که مدل های‌اسپلیت برای هر روز آلوده سه مسیر طوفان‎ 
گردوغباری را پیش‌بینی می‌کند و برآیند مسیرهای پیش‌بینی شده به‌عنوان کانون مولد گردوغبار معرفی می‌شوند.‎ 

الگوی داخل اسستان: با توجه به برآیند تحلیل‌های انجام شده جهت‌های وزش باه سطح اثرگذاری رخدادهای 
گردوغبار و مدل‌سازی انجام شده موارد معدودی از رخدادهای گردوغباری استان همدان دارای منشأ داخل استانی 


هستند. عمده‌ترین احتمال وقوع الگوی داخل استان در شکل (۳) نمایش داده شده است. 


ms 
3 60 Kilometers 


شکل ۳- الگوی حرکت ذرات در مدل های‌اسپلیت برای روز وقوع طوفان گردوغبار در داخل استان (الگوی 
اول) 


می‌گیرند. مدل‌سازی های‌اسپلیت نشان داد که بخش اعظم گردوغبارها از استان‌های واقع در غرب و جنوب غربی 


سال یازدهم مسیریابی حرکت گردوغبار بر اساس غلظت .... Mo‏ 


استان همدان وارد می‌شوند. بخش دیگری نیز از استان‌های واقع در شرق و شمال شرقی محدوده مطالعاتی و نواحی 


مرکزی ایران sa‏ می‌گیرند. الگوی LS‏ این حرکت‌ها در شکل (E)‏ نشان داده شده است. 
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شکل £- الگوی حرکت ذرات در مدل های‌اسپلیت برای روز وقوع طوفان گردوغبار در خارج از استان (الگوی 


دوم) 


الگوی خارج از کشسور: در چند سال اخیر به دلیل وفور احداث سد در کشورهای ترکیه سوریه و اران و همچنین 
تغییرات رژیم آبی و همچنین خشکسالی حاکم در منطقه سبب خحشک شدن بسیاری از تالاب‌هاء از بین رفتن پوشش 
گیاهی و خشک شدن منطقه گردیده است. از طرفی در اثر تماس باد با خحاک» گردوغبار ب‌صورت ریزدانه از سطح 
NS TL‏ ل کی رکد پررسی سیر هان نخر کت تراک کرو غاز با نتفای از pilag‏ 
فان واه yea‏ سین Le o T‏ ارج Gal odds‏ مان ا 2 oaa A MIS‏ رل ارب 
که beast‏ از غراق اران کل a‏ کل Lie calls al aA‏ هراق تردن با cala dou‏ اا د &4 
باشد که از طریق مرز عراق وارد ایران و ble‏ وارد استان همدان می‌شوند. مسیر دوم از جنوب غرب گردوغبار را از 
شبه‌جزیره عربستان و يا بخش‌های جنوبی کشور عراق (در مجاورت استان خوزستان) منتقل می‌کنند و دلیل اصلی 
این طوفان‌ها می‌تواند خشک شدن تالاب‌های مرزی و عدم مدیریت صحیح مناطق خشک در کشورهای همسایه 


unl 


141 جغرافیا و مخاطرات محیطی شمارة اول 


راهنما 
ATH |‏ — 
Y ۷‏ - | 
sem ۴‏ 
Cote!‏ 
انان مدان که Kilometérs‏ 0 


شکل ۵- الگوی حرکت ذرات در مدل های‌اسپلیت برای روز وقوع طوفان گردوغبار در خارج از ایران (الگوی 


سوم) 


-٤‏ نتیجه گیری 

به‌طورکلی س‌الانه استان همدان در معرض طوفان‌های متعددی با Lites‏ محلی منطقه‌ای و همچنین فرامنطقه‌ای 
است. وقوع طوفان در مناطق خشک ازجمله استان همدان سبب افزایش ذرات معلق نسبت به حد آستانه شده و این 
موضوع اثرهای منفی بر سامت انسان دارد؛ بنابراین» یکی از راهکارهای اساسی در کنترل آلودگی‌هایی نظیر 
طوفان‌های گردوغباری شناسایی منبع آلودگی است. ازجمله روش‌های شناسایی برخحاستگاه‌ها استفاده از روش‌های 
سنجش ازدور و مدل‌های عددی مسیر انتقال جریان هوا است. سنجش ازدور با فراهم آوردن تصاویر با پوشش وسیع 
و تکرار زیاد مطالعات گردوغبار را با سرعت و دقت زیاد عملی نموده است و در حال حاضر مناسب‌ترین ابزار برای 
بررسی این پدیده است؛ به‌طوری‌که اطلاعات زیادی را در زمینه مشخصات طیفی گردوغبار ازجمله شناسایی 
گردوغبان پایش و ردیابی مسیر گردوغبارء شناسایی کانون‌های ایجاد گردوغبار و برآورد شدت گردوغبار فراهم 
UR a‏ در این پژوهش به استفاده از سنجش ازدور و اطلاعات هواشناسی جهت مسیریایی طوفان‌های گردوغباری 
غرب کشور پردانخه شید در dag‏ اضر از sluts pAOD yati‏ کانرن‌های گردوغبار teal‏ شد که 
eal in cola ol‏ آن را فان cnl‏ سیر و راد OA‏ و SOW DARE Jas‏ 
UE TERT AST a SEE RE LUTEA‏ 
معلق ب وف کر غار bo‏ دیکن قاویایی dials oos Reng‏ که LAOD cleats‏ 


استفاده از روش مستقیم» توزیع مکانی بسیار خوب الگوی آئروسل محلی در مقایسه با روش‌های دیگر را دال 


tas,‏ اما مشاه توهش عیاش و کی E‏ سب py‏ ماد AUS Ee fsa:‏ ازسمله وان فکیک پان 
ممکن است در مقایسه با تصاویری با توان تفکیک Ub‏ در شناسایی دقیق کانون تولید ریزگرد کارآمد نباشد. همچنین» 
در این پژوهش مشابه گزارش ناسا مقدار ۰/۳ به‌عنوان آستانه گردوغبار در نظر گرفته شده است؛ درصورتی که در 
پزوهش دانیالی و همکاران (۱۳۹۷) عمق اپتیکی بالاتر از ۰/۵ به‌عنوان معیاری از وقوع طوفان گردوغباری در نظر گرفته 
شده است. درمجموع ۲۵۰ رخداد گردوغبار در بازه زمانی موردمطالعه بررسی شد که این رخدادها از موارد وسیع که 
تمام یا بخش‌های عمده‌ای از استان را در برمی‌گرفت تا مواردی با دامنه بسیار محدود را شامل شده است. نتایج به 
دست آمده در این پژوهش نیز با نتایج سایر پژوهش‌های میری و همکاران YAT)‏ دانیالی و همکاران (۱۳۹۷) که 
سه basal aso ec‏ سل سانش cl ca cles‏ کر eot alus‏ دارو و نخان تلم این 
موضوع است که این مدل در مسیریابی گردوغبار موفق عمل کرده است؛ ازاین‌رو نتایج حاصل از حروجی مدل‌سازی 
های‌اسپلیت نیز OLY‏ داد که منشأً گردوغبارهای استان همدان در طوفان‌های بررسی شده دارای چهار مسیر کلی 
است. عمله‌ترین بخش ورودی استان» از محور جنوب و جنوب غرب استان است که منتقل کنندة گردوغبارهای با 
منشأً حارج از استان و خارج از کشور است. در این پژوهش نیز با توجه به تحقیقات مشابهی ازجمله فلاح زوزولی و 
ENTIS S ur. T‏ کم رة اردق ر ران کان هان اتی 
ورود گردوغبار به داخل کشور هستند که متعاقباً استان‌های مرزی متأثر از این پدیده می‌گردند. دیگر مسیر ورودی از 
کر gd‏ کا بد رل فر وان اسان وو باعل ان ورد س شیف ورک عمد وهای 
گردوغبار Lites‏ خارج از استان دارند dy‏ می‌توان مواردی (حداکثر ۱۸ مورد) در داخل محدوده مطالعاتی دانست. 
IS sha‏ بخش قابل توجهی از کانون‌های شناسایی‌شده در خارج از مرز سیاسی استان و در داخل بافر ۱۵ کیلومتری 
در استان‌های مجاور قرار گرفته است که به‌هرحال تأثیر آن در بروز پدیده گردوغبار در داخل استان همدان انکارناپذیر 
است. ازاین‌رو بخش‌هایی از استان‌های کردستان» مرکزی خوزستان و همچنین استان زنجان از مهم‌ترین بخش‌های 
US SE‏ در بروز و انتقال گردوغبار به داخل استان هستند. در داخل استان نیز بخش‌هایی از شسهرستان‌های فامنین» 
ملایر و رزن نیازمند مدیریت جدی جهت JES‏ پدیده گردوغبار در آینده هستند. مشاهدات میدانی نیز نشان از 
lesa‏ اا ی وهای اقاس هو غات د غا "TRU UL E TII‏ 
اماق هرا اکن و cales‏ نکی از coe ela o‏ در مطالعات wate‏ ر فان‌های eol Jess;‏ و از 
s aus‏ مها سل ad eoe‏ و قط فرت cal‏ فرش شاه اد ji gallo pa dbi‏ 
E TE slik‏ 
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